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Summary 

The  P  system  is  a  general  distributed  model,  highly  parallel  and  based  on  the  notion  of  a 
Membrane  structure.  This  investigation  addressed  the  development  of  a  model  for  P 
systems  with  Active  Membranes  using  the  Membrane  Creation  technique.  The  usefulness 
of  the  model  was  shown  by  applying  the  model  to  solve  the  Hamiltonian  Path  Problem 
(HPP)  for  Undirected  Graphs.  The  relevant  algorithm  for  the  problem  is  presented. 


Introduction 

P  systems  are  a  class  of  distributed  parallel  computing  models  introduced  in  (1),  inspired 
by  the  way  live  cells  process  chemical  compounds,  energy,  and  information.  In  short,  in 
the  regions  delimited  by  a  membrane  structure,  one  places  multi-sets  of  objects,  which 
evolve  according  to  evolution  rules  associated  with  the  regions.  A  computation 
consists  of  transitions  among  system  configurations.  The  result  of  a  halting  computation 
is  the  vector  of  the  multiplicities  of  objects  present  in  the  final  configuration  in  a 
specified  output  membrane  or  of  objects  which  leave  the  external  membrane  of  the 
system  (the  skin  membrane)  during  a  computation. 

The  P  systems  can  be  of  different  forms,  such  as: 

1)  P  system  with  Labeled  Membranes 

2)  P  system  with  Polarized  Membranes 

3)  P  systems  with  Active  membranes. 

In  this  investigation  we  addressed  P  systems  with  Active  membrane  defined  in  (2,  3). The 
P  systems  with  Active  Membranes  can  be  classified  in  two  fonns: 

1)  Passive  P  systems,  and 

2)  Active  P  systems. 

The  P  systems  addressed  in  (2,  3)  are  Passive  systems.  Here  the  membrane  structure 
consists  of  m  different  membranes  and  during  the  computation  this  number  may  either 
decrease  (by  dissolving  membranes)  or  increase  (by  dividing  the  existing  membranes)  but 
this  number  (  with  respect  to  the  labels  of  membranes)  is  always  less  than  or  equal  to  m. 
The  degree  of  a  passive  P  system  is  the  number  of  initial  membranes  present  in  the 
system. 

Sometimes  it  may  happen  that  different  types  of  membranes  may  be  increased  in  a 
membrane  system  by  creating  new  membranes  whose  labels  are  different  from  the 
existing  ones.  Such  systems  are  called  Active  P  systems. 

Membranes  are  created  continuously  in  biology,  for  instance  in  the  process  of  Vesicle 
mediated  transport.  Because  one  of  the  roles  of  membranes  is  to  keep  the  molecules  of 
a  compartment  close  to  each  other  in  order  to  facilitate  their  reactions,  when  a 
compartment  becomes  too  large,  it  often  happens  that  new  membranes  appear  inside  it, 
more  or  less  spontaneously  or  during  biological  evolution.  The  process  of  creating  new 
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membranes  is  also  known  as  “Autopoiesis”  and  the  details  of  the  process  are  addressed 
in  (4,  5). This  investigation  considered  the  creation  of  membranes  under  the  influence  of 
existing  objects.  We  abstracted  the  operation,  using  it  in  an  idealized  way,  mainly  as  a 
tool  for  better  structuring  the  objects  of  a  given  membrane  system. 

Membrane  creation  can  be  introduced  both  for  systems  with  Symbol  objects  (6-9)  and 
for  systems  with  String-objects  (10).  We  considered  the  case  of  symbol  objects. 

The  P  system  with  Active  Membranes  was  investigated  in  the  2003  Visiting  Faculty 
Research  Program  (VFRP)  research  Project  (11).  The  model  was  enhanced  in  the  2004 
VFRP  research  project  and  the  algorithms  for  the  Satisfy-ability  problem  (SAT)  and  the 
HPP  problems  are  given  in  the  Appendix.  In  this  project  the  Membrane  Creation 
Techniques  were  applied  to  develop  the  P  System  model. 

The  following  sections  describe  the  model  and  its  application. 


The  Model 

A  P  system  with  Membrane  Creation  of  degree  (m,  n),  n  >  m  >  1,  is  a  construct: 


n  =  (V,  T,  C,  p,  w0,  w,, . ,  W(m-i),  Ro,  Ri, . ,  R(n-i) ), 

Where: 

1.  V  is  an  alphabet;  it  consists  of  both  productive  and  non-productive  objects; 

2.  T  c  V,  is  the  output  alphabet; 

3.  C  n  V  ^  cp,  is  the  set  of  catalysts; 

4.  p  is  a  membrane  structure  consisting  of  m  membranes,  with  the  membranes  and 

the  regions  labeled  in  a  one-  to-one  manner  with  elements  in  a  given  set;  here  the 
labels  are  used  as  0  (for  the  skin  membrane),  1, . ,  (m-1); 

5.  Wi ,  0  <  i  <  (m-1),  are  multi-sets  of  objects  over  V  associated  with  the  regions 
0,  1,....,  (m-1)  of  p; 

6.  Ri5,  0  <  i  <  (n-1),  are  finite  set  of  evolution  rules  over  V. 

An  evolution  rule  is  of  two  types: 

(a)  If  a  is  a  single  non-productive  object,  then  the  evolution  rule  is  in  the  form  a  — >  v 
or  c  a  — >  c  v,  where  ceC,  a  e  (V  -  C),  v  =  v  or  v  =  v  8  or  v  =  v  x,  where  v  is  a 
multi-set  of  objects  over  ((V-C)  x  (  here,  out))  u  ((V-C)  x  ( in,  I  1  <  j  <  (n-1))), 
and  8,  t  are  special  symbols  not  in  V. 
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(b)  If  a  is  a  single  productive  object,  then  the  evolution  rule  is  in  the  form 

a  — >  v  ] ,.  or  c  a  — >  c  [,.  v  ]  . ,  where  ceC,  a  e  (V  -  C),  v  is  a  multi-set  of 

objects  over  (V-C).  The  rule  of  the  form  a  — >  [(.  v  ]  means  that  the  object 
identified  by  a  is  transformed  into  the  object  identified  by  v,  surrounded  by  a 
new  membrane  having  the  label  i.  No  rule  of  the  form  a  — >  [0  v  ]  0  can  appear  in 

any  set  R, .  During  a  computation  the  number  of  membranes  can  increase  or 
decrease. 

The  membrane  structure  and  the  multi-sets  in  II  constitute  the  initial  configuration  of  the 
system.  One  can  pass  from  a  configuration  to  another  one  by  using  the  evolution  rules. 
This  is  done  in  parallel:  all  objects,  from  all  membranes,  which  can  be  the  subject  of  local 
evolution  rules,  should  evolve  simultaneously.  A  rule  can  be  used  only  if  there  are  objects 
which  are  free  at  the  moment  one  checks  its  applicability. 

The  application  of  a  rule  c  a  — >  c  v  in  a  region  containing  a  multi-set  w  means  to 
remove  a  copy  of  the  object  a  in  the  presence  of  c  (catalyst),  provided  that  such  copies 
exist,  then  follow  the  prescription  given  by  v:  If  an  object  appears  in  v  in  the  form  (a, 
here),  then  it  remains  in  the  same  region;  if  it  appears  in  the  form  (a,  out),  then  a  copy  of 
the  object  a  will  be  introduced  in  the  region  of  the  membrane  placed  outside  the  region 
of  the  rule  ca-^cv  ;  if  it  appears  in  the  form  (a,  in; ),  then  a  copy  of  a  is  introduced  in 
the  membrane  with  index  i,  if  such  a  membrane  exists  inside  the  current  membrane, 
otherwise  the  rule  can  not  be  applied.  If  the  special  symbol  x  appears,  then  the  thickness 
of  the  membrane  which  delimits  the  region  in  question  is  increased  by  1 .  Initially,  all 
membranes  have  the  thickness  1 .  If  a  rule  in  a  membrane  of  thickness  1  introduces  the 
symbol  x,  then  the  thickness  of  the  membrane  becomes  2.  A  membrane  of  thickness  2 
does  not  become  thicker  by  using  further  rules  which  introduce  the  symbol  x,  but  no 
object  can  enter  or  exit  it.  If  a  rule  which  introduces  the  symbol  5  is  used  in  a  membrane 
(including  the  skin  membrane)  of  thickness  1,  then  the  membrane  is  dissolved.  If  the 
membrane  had  thickness  2,  then  it  returns  to  thickness  1 .  When  ever  the  skin  membrane 
is  dissolved,  the  whole  membrane  system  will  be  destroyed.  If  at  the  same  step  one  uses 
rules  which  introduce  both  5  and  x  in  the  same  membrane,  then  the  membrane  does  not 
change  its  thickness.  No  object  can  be  communicated  through  a  membrane  of  thickness 
two,  hence  rules  which  introduce  commands  out  and  in,  requesting  such  communications, 
can  not  be  used.  The  communication  has  priority  over  changing  the  thickness.  If  at  the 
same  step  an  object  should  be  communicated  and  a  rule  introduces  the  action  x,  then  the 
object  is  communicated  and  “after  that”  the  membrane  changes  the  thickness.  When 
applying  a  rule  a  — >  v  ] ,  in  a  region  j,  a  copy  of  a  is  removed  and  a  membrane  with 

the  label  i  is  created,  containing  the  multi-set  v,  inside  the  region  of  membrane  j.  We  can 
never  create  a  skin  membrane.  Similarly,  when  applying  a  rule  c  a  — >  c  [,.  v  ] ,  in  a 

region  j,  a  copy  of  a,  in  presence  of  c  (catalyst),  is  removed  and  a  membrane  with  the 
label  i  is  created,  containing  the  multi  set  v,  inside  the  region  of  membrane  j. 

A  sequence  of  transitions  between  configurations  of  a  given  P  system  II  is  called  a 
Computation  with  respect  to  IT.  A  computation  is  successful  if  and  only  if  it  halts,  that 
is,  there  is  no  rule  applicable  to  the  objects  present  in  the  last  configuration.  The  result  of 
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a  successful  computation  is  \|/T(w),  is  denoted  by  Ps(Il)  (from  “Parikh  set”)  and  we  say 
that  it  is  generated  by  II.  (Note  that  we  take  into  account  only  the  objects  from  T). 

The  model  was  applied  to  solve  the  following  problem: 


The  Hamiltonian  Path  Problem  (HPP)  for  Undirected  Graphs. 

Let  G  =  (U,  E)  be  an  undirected  graph  with  n  nodes  (n  >  2),  where  U  is  the  set  of  nodes 
and  E,  the  set  of  edges. 

The  Hamiltonian  path  problem  determines  whether  or  not  there  exists  a  Hamiltonian  path 
in  G,  that  is,  whether  or  not  there  exists  a  path  that  passes  through  all  the  nodes  in  U 
exactly  once.  The  HPP  for  undirected  graphs  is  known  to  be  Non  Deterministic 
Polynomial,  NP-Complete. 

Let  U  =  (ab  a2, . ,  an ).  We  now  construct  a  P  System  with  membrane  creation  of  degree 

(1,  n  +  1),  where  1  is  the  initial  number  of  membranes  &  (n+1)  is  the  total  number  of 
membranes  used  in  the  computation. 


n  =  (V,  T,  C,  p,  Wo,  R0,  Ri, . ,  Rn ), 

Where  V  =  (  a, ,  a,  ,  f ,  f  ,  d, ,  d,  ,  c,J  I  1  <  i  <  n,  0  <  j  <  n  )  u  (  Y  )  u 


(  h  I  0  <  I  <  3n  );  (  here  a;  are  productive  objects  ). 


T  =  (Y); 

C  =  cp; 

H  [o  ]  o ’ 


Wo  —  ( t0 ,  ai  ,  a2 , . ,  an ); 

The  set  R0  contains  the  following  rules: 

1 .  f  ->  ti+i ,  0  <  i  <  3n  -  1 ; 

2.  t3n  — » 5; 

3.  c/  — *■ X  ,  1  <j  <n  -  1,  V,.; 

4.  c"  — >  ( Y,  out ),  V. ; 

5.  ai->  [,.  f  ,  Cj° ,  ajl? . ,ajk  ].,  1  <  i  <  n; 

(  ai  will  create  a  new  membrane  with  label  i  and  transform  into  a  multi-set  of  objects 
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f  ,  Ci° ,  aji, . ,  ajk  ,  where  ji , . ,  jk  are  vertices  adjacent  to  vertex  i  ). 

The  set  R; ,  1  <  i  <  n  ,  contains  the  following  rules: 

I-  aj -*■  [  ^ dj ,a j!, . , ajk  ]j,  1  <  j  <  n; 

(  a j  will  create  a  new  membrane  with  label  j  and  transform  into  a  multi-set  of  objects 
fj  ,  dj  ,a  ji, . ,  ajk  ,  where  j, , . ,  jk  are  vertices  adjacent  to  vertex  j ). 

2.  dj-di; 

3.  f;  -» (  f ,  inj,  )  (  f ,  inj2 )  . (  f  ,  injk ); 

(Here  ji , . ,  jk  are  vertices  adjacent  to  vertex  i  ). 

4.  fj  ->  (  fj ,  injx  )  (  fj ,  inj2 )  . (  fj ,  injk )  ,  Vy.  *  i ; 

(Here  j] , . ,  jk  are  vertices  adjacent  to  vertex  i  ). 

5.  a  j  — >  aj ; 

6.  d;  — >  C°i  ; 

7.  f— >Xt; 

8.  ckj  -»  (  Cjk+1 ,  out ) ,  k  >  0  ; 

N,  lin+D(  n  )  =  (  Ym  ,  if  the  given  graph  has  m  Hamiltonian  paths.  ) 

=  cp  otherwise. 

The  system  works  as  follows: 

Initially  the  skin  membrane  contains  all  the  productive  objects  ax  ,  a2 , . ,  an,  where  each 

productive  object  a; ,  by  using  the  rule  a,  — >  [;  f  ,  c,° ,  aJb . ,  ajk  ],.,  will  create  a 

membrane  with  index  i  and  transform  into  the  multi-set  of  objects 

f  ?  Cj  ?  4j  1’ . '  ajk. 

In  membrane  i,  a  productive  object  ajwill  create  a  membrane  with  index  j  and  transform 

into  a  multi-set  of  objects  fj'',  d,  ,aiu . ,  a  jk.  In  membrane  i,  f  j  will  be  replaced  with  fj , 

d  j  will  be  replaced  with  d, ,  and  a  j  with  aj .  In  the  next  step,  d,  will  be  replaced  with  c°j ,  aj 
will  create  a  membrane  by  using  the  rule  a,  — ►  [(  fj”,  d,  ,a jX, . ,  ajk  ]  and  the  process 

defined  above  will  be  repeated.  When  ever  a  new  membrane  is  created,  all  fj’s  and  f  fs, 
from  its  parent  membrane,  will  move  into  that  membrane.  The  role  of  fj  is  to  increase  the 
thickness  of  the  membrane  j,  so  that  no  object  will  be  sent  out  from  that  membrane.  This 
will  eliminate  the  effect  of  the  path,  which  contains  multiple  copies  of  the  same  node,  on 
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the  result.  For  every  rewriting  step,  Cjk  will  be  rewritten  as  Cjk+1  and  sent  out  of  the 
membrane  i  provided  that  the  thickness  of  the  membrane  i  is  1 .  In  the  skin  membrane,  c", 
(indicates  that  there  exists  a  tour  of  path  length  n)  will  be  replaced  with  Y  and  sent  out.  A 
counter  f  is  used  in  the  skin  membrane,  so  that  for  every  rewriting  step  it  will  be 
incremented.  When  ever  the  counter  reaches  t3n  ,  indicating  the  end  of  the  computation, 
the  counter  will  dissolve  the  whole  membrane  system  by  dissolving  the  skin  membrane. 
If  there  exists  a  tour  of  length  n,  by  traversing  all  nodes  in  the  graph,  then  we  get  Y  from 
the  system. 

Time  Complexity:  This  algorithm  takes  3n+l  steps  for  generating  the  output. 


The  Algorithm:  Hamiltonian  Path  Problem  for  Undirected  Graphs. 

Let  us  consider  the  graph  G  =  (U,  E),  with  U  =  {a,,  a2,  a3}  as  shown  Fig  1  below: 


Figure  1:  The  Undirected  Graph  (HPP) 

We  now  construct  a  P  system  with  membrane  creation  of  degree  (1,  4)  as: 

II  =  (V,  T,  C,  p,  Wo,  R0,  . ,  R3 ), 

Where  V  =  (  at,  a,  ,  f,  f  ,  d, ,  d,  ,  c,J  I  1  <  i  <3,  0  <  j  <  3  )  u  (f  I  0  <  i  <  9)  u  (Y); 
T  =  (Y); 

C  =  cp; 

h  _  [o  ]  0  ’ 
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Wo  —  ( t0,  ai ,  . ,  a3 ); 

and  R; ,  0  <  i  <3, 

1 .  [  o  t0  a  i  a2  a3  J  0 

2.  [0  ti  [j  f, Cl°a2  ]1[2f2'c2°a1a3  ]  2  [3  f3  c3  0  a2  ]  3  ]0 

3.  [0 12  cd  c2‘  c3‘  [j  fj' '  [2  f2  d2  a,  a3  ]2  ]  i  [2  f2  [x  f/  d/  a2  ]  3  [3  f3  d3  a2  ]3  ]2 

[3  f3  [2  f2  d2  a!  a3  ]  2  ]  3  ]  0 

4.  [0  t3  [1[2f1f2’d2a1a3  ]2  ],[2[1fi  f2d,a2  ] ,  [3  f2  f,  d3  a2  ]3  ]2 
[3[2  f 2  ^3  d2  aj  a3  ]  2  ]3  ]0 

5-  [0  [ i [ 2  fi  f2  c2°  [,  fi  d,  a2  ]J3  f3  d3  a2  ]3  ]2  ] , 

[2[i  fi  f2C!°[2  f2  d2  a/  a3  ]2  ] ,  [3  f2  f3  c3°[2  f2  d2  a!  a3  ]2  ]3  ]2 

[3[2  f2  f3C2°[1  fi  di  a2  ]j[3f3  d3  a2  ]3  ]2  ]3  ]0 

6-  [o  tsti  c2‘  [2[x  fi  f2f|  d,  a,  ] ,  [3  f,  f2  f3  d3  a,  ]3  ]2  ], 

[2  Ci1  c3'  [t  [2  f,  f2  f2  d2  ai  a3  ]  2  ]{ [3  [2  f2  f3  f2  d2  a;  a3  ]  2  ]  3  ]  2 

[3c2‘[2[i  f2  f3  f/  dj  a2  ]j  [3f2f3f3  d3a2  ]3  ]2  ]3  ]0 

7.  [0  t6  c,2  c22c22  c32  [,  [2  (j  f2  f/  Ci°  [2  f2  d2  a/  a3'  ]  2 ), 

[3  f,  f2  f3  c3°[2  f2  d2  a,  a3  ]2  ]3  ]2  ], 

[2[i  (2  fi  f2  c2°  [j  fi  di  a2  ]  j  [3  f3  d3  a2  ]  3)2  ]  l  [3  (2  f3  f2  c2°  [3  fi  di  a2  ]  3 
[3  f3  d3  a2  ]  3)2  ]  3  ]  2 

[3 [2 [1  f2  f3  fi'  Ci°[2  f2  d2  a!  a3'  ]2  \  (3  f2  f3' c3°  [2  f2  d2  aj  a3'  ]2  )3  ]2  ]3  ]0 

8.  [0  t7  [,  [2  c3‘  (1  Cl°  [2  f2  f,  f2  d2  a,  a3  ]2  )x  [3[2  f,  f2  f,  f2  d2  a,  a3  ]2  ]3  ]2  \ 
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[2  [i  (2  c2°  [j  fif2fi  d,a2  ]j[3  f;  f2f3  d3  a2  ]3  )2  1[3(2  c2°  [,  f3  f2  f,'  d,  a2  ], 

[3  f3  f2  f3  d3  a2  ]  3  )2  ]3  ]2 

[3  [2  ci‘  [1  [2  f2  f?  fi  f2  d2  aj  a3  ]2  \  (3c3°  [2  f2  f3  f2  d2  a;  a3  ]2)3  ]2  ]3  ]0 

9.  [0  tg  [{  c32  [2 3(  c !°  (2  f,  f2  c2°  [t  f,  d/  a2  ]j[3  f3  d3  a2'  ]3  )2  )3  [3  (2  f,  f3  f2  c2° 

[1  fi  di  a2  ]  j 

[3  fs  d3  a2  ]3  )2  ]3  ]2  ]j 

[2[i  (2  c2°  (1  f2  fi  Ci°  [2  f2  d2  aj  a3  ]2)i  [3  fi  f2  f3  c3°  [2  f2  d2  ai  a3  ]2  ]3  )2  ]t 
[3  (2  c2°  [j  f3  f2  fi  ci°  [2  f2  d2  a)  a3  ]2  ]:  (3  f2  f3  c3°  [2  f2  d2  a]  a3  ]2  )3  )2  ]3  ]2 

[3  c,2  [2[t  (2  f3  f2  c2°  f,  d/  a2  ]j[3  f3  d3  a2  ]3  )2  ],  (3  c3°  (2  f3  f2  c2°  [t  f/  d/  a2  ]; 

[3  f3  d3  a2  ]3  )2  )3  ]2  ]3  ]0 

10.  [0 19  C]3  c33[j[2  (iCi°(2c2°  [,  f,  f2  f,  d,  a2  ],[3  fjf2f3  d3a2  ]3)2), 

[3  (2  c2°  [j  f,  f3  f2  f,  d,  a2  \ 

[3  f  1  f3  f2  f3  d3  a2  ]3  )2  ]3  ]2  ]j 

[2[i  (2  c2°  c3‘  (1  Ci°  [2  f2  fi  f2  d2  &i  a3  ]2 ) \  [3[2  f!  f2  f3  f2  d2  a;  a3  ]2  ]3  )2  \ 

[3  (2  c2°  Ci'  [j[2  f3  f2  fi  f2  d2  aj  a3  ]2  ]j  (3  c3°  [2  f2  f3  f2  d2  ai  a3  ]2  )3  )2  ]3  ]2 

[3  [2  [1(2  c2°  [j  f3  f2  f,  d,  a2  ]j[3  f3f2f3  d3a2  ]3)2  \  (3  c3°  (2  c2°  [,  f3  f2  f,  d,  a2  ]j 
[3  f3  f2  f3  d3  a2  ]3  )2  )3  ]2  ]3  ]0 

11.  Y2 

2 

Y  indicates  that  two  (2)  Hamiltonian  paths  exist  in  the  given  graph. 

The  paths  are:  aj  -  a2  -  a3  and  a3  -  a2  -  a!. 

Note:  Here  a  different  notation  (  )  is  used  for  membranes  of  thickness  2. 
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Conclusions 

The  mathematical  model  for  the  P  system  with  Active  Membranes  was  developed.  The 
model  was  applied  to  solve  the  Hamiltonian  Path  Problem  for  Undirected  Graphs  and  the 
relevant  algorithms  developed. 

Recommendations  for  Follow  on  Research 

The  following  Research  Topics  are  identified  for  further  investigation: 

1)  Development  of  a  Simulation  model  in  Digital  Computer  for  P  Systems 
with  Active  Membrane, 

2)  Design  and  develop  experimental  techniques  to  implement  the  model  in 
Bio-Chemical  Media. 
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Appendix  A:  Algorithm  to  solve  SAT  problem 

The  computational  steps  for  the  construction  of  a  P  system  for  the  propositional  formula: 

P  =  (xi  V  x2  V  x3)  A  (~X|  V  ~x2  V  x3)  A  (xi  V  ~x2  V  ~x3)  with  n=3  and  m=3 
are  shown  below: 

Step  0: 

The  initial  configuration  of  the  system  is  given  as: 

[4  [3  [2  [1  [0  coai  a2a3  ]  °  ]  1  ]  2  ]  3  ]  ° 

Step  1: 

[4  [3  [2  [1  [0  citia2a3  ]  0+  [0  c1f1a2a3  ]  -  ]  °  ]  °  ]  °  ]°4 


;  Lofi  Jo 

Step  2: 

Membrane  1  is  divided  due  to  the  presence  of  the  two  copies  of  membrane  0  with 
opposite  polarizations  using  the  rule: 

[,  [»  E  [.  t  If  -  [,  [.  12  If  [,  [.  If  f 

The  counter  ci  is  replaced  by  c2.  No  new  truth  value  is  introduced  at  this  step  since 
membrane  0  is  not  of  neutral  polarity.  We  get, 

[4  [3  [2  [1  [0  C2  fi  a2  a3  Jo  ] [j  [0  c2  fi  a2  a3  Jo  ] ,  ]  2  ]  3  J  4 


The  rules  applied  in  parallel  are: 

[0  Co  ci  ]  °  and 

the  division  rule,  [0  a3  ]  °  — >  [0 1, 
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Step  3: 

Membrane  0  is  divided  and  the  truth  assignment  of  the  variable  a2  also  takes  place 
according  to  the  following  rule: 

\oai  ]o  *  [oC  ]o [0 f 2  ]o 

The  counter  C2  is  replaced  by  C3.  Membrane  2  is  also  divided  due  to  the  presence  of  two 
membranes  1  with  opposite  polarity  by  the  rule: 

[J.  id,  t  E  -  u,  is  e  u,  is  t 

We  get: 

[4  [3  [2  [1  [0  c3  h  b  &3  ]  0  [0  c3  ti  £2  a3  ]  0  ] ,  ]2 

[2  [1  [0  c3  ft  k  a3  ]  0  [0  C3  fi  f2  a3  ]  o  ]  1  ]2  ]  “  ]  4 


Step4: 

Membrane  1  is  divided  since  it  contains  two  membranes  0  of  opposite  polarity  using  the 
rule: 


[,  [„  l  [.  Is  ]?  -  [,  [„  IS  k  [,  [„  IS  Ir 

The  counter  c3  is  replaced  by  C4.  Membrane  3  would  also  be  divided  since  it  contains 
membranes  2  of  opposite  polarity  using  the  rule: 

f  f  i+  r  1  i°  r  r  i°  i°  [  [  i°  i° 

L  m  Lm—\  J  m— 1  L  m—l  J  m— 1  Am  L  m  Lm-1  J  m— 1  Am  L  m  L  m-1  J  m- 1  J  m 


since  here  m=3.  Thus,  we  get, 


[4 


[3  [2  [1  [0  c4  ti  b  a3  is  ]:[,[  0  C4  ti  f2  a3  ]  q  ]  j 

[j  [2  [,  [0  <=4  f.  t2  a3  ];  ];  [,  [,,  c4f,  t;  a.  ];; 


IS  1 


]r  IS  IS  ] 


0 

4 


Step  5: 

Membrane  0  is  divided  again  and  there  occurs  truth  assignment  of  the  variable  a3.  Also, 
membrane  2  is  divided  due  to  the  presence  of  two  opposite  polarity  membrane  1 .  Counter 
C4  is  replaced  by  C5.  Thus,  we  get: 
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L  [)  [j  [|  [»  C5  ti  t2  t,  ]„  [„  C5  ti  t2  fj  ]„  ]  f  ]:  [,  [,  [„  C5  t,  f2  t3  1,  [0  C5  t,  f2  f j  ],, 

]f  t  IS 

[;  [2  [,  [« cs  fi  t,  t3  ];  [„  c5  f,  t,  fj  t  ]  ;■  ]:  [,  [,  [„  cs  f,  t2 t3  ],  [0  c.  f,  f2  f3  ]„ 

]f  t  IS  IS 

Step  6: 

Membrane  1  is  divided  by  the  rule 

[,  [»  E  [»  ]„  IS  -  [,  [»  IS  ];  [,  [„  12  ]r 

Membrane  3  is  divided  by  the  rule 

[  [  i+  r  l  i°  _h.  r  r  i°  i°  r  [  1°  1° 

L  m  Lm-1  J  m— 1  Lm-1  J  m- 1  Am  L  m  L  m-1  J  m— 1  Am  L  m  L  m-1  J  m— 1  J  m 

since  here  m=3.  The  counter  C5  is  replaced  by  C6.  Thus  we  get, 

[4  [3  [2  [1  [0  c6  ti  t2 13  ]  0  ]  1  [1  [0  c6  ti  t2  f3  ]  0  ]  i  ]  2  ]  3 

[3  [2  [1  [0  c6  ti  f2 13  ]  0  ]  t  [1  [0  c6  ti  f2  f3  ]  [J  ]  T  ]  2  ]  3 

[3  [2  [l  [0  C6  fl  t2  t3  ]  0  ]  1  [j  [0  C6  fl  t2  f3  ]  0  ]  I  ]  2  ]  3 

[3  [2  [l  [0  C6  fl  f2  t3  K  ]r  [,  [0C6f1f2f3  ]°  ]~,  ]“  ]3°  ]“ 

Step  7 : 

No  more  truth  assignments  are  possible,  so  membrane  0  cannot  divide  any  further.  The 

counter  C6  is  replaced  by  C7.  Membrane  2  is  divided  and  we  get, 

[4  [3  [2  [l  [0  c7  tl  t2  t3  ]  0  ]  1  ]+2  [2  [l  [0  C7  tl  t2  f3  ]  0  ]  l  l  ]  3 

[3  [2  [l  [0  C7  tl  f2  t3  ]  0  ]  1  ]+2  [2  [l  [0  C7  tl  f)  f?  ]  0  ]  I  ]2  ]  3 

[3  [2  [l  [0  C7  fl  t2  t3  ]  0  ]  1  ]+2  [2  [l  [0  C7  fl  t2  f?  ]  0  ]  1  L  ]  3 

[3  [2  [l  [0  C7  fl  f2  t3  ]°  ]?  I  [2  [,  [0  C7  fl  f2  f3  ]S  ]?  ]2  ]“  ]l 
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Step8: 

Membrane  3  is  divided  by  the  rule, 

[  f  i+  r  l  i°  _>  r  r  i°  i°  r  r  1°  i° 

L  m  Lm-1  J  m— 1  L  m—l  J  m— 1  J  m  L  m  L  m—l  J  m- 1  J  m  L  m  L  m-1  J  m—l  J  m 

The  counter  C7  is  replaced  by  c8.  We  get, 

[4  [}  [2  [l  [0  c8  ll  *2  t3  ]  0  ]  1  ]  2  ]  3  [3  [2  [l  [0  C«  tl  t2  f3  ]  0  ]  1  ]  2  ]  3 

L  L  [,  [„ oB t, f2 13  ];  ]?  ];  ];  [,  [2 [,  [„ c8 1, f2 f,  ]2  ]f  ]2  ]; 

[3  [2  [l  [0  C8  fl  U  t3  ][!  ]  2  ]2  ]2  [3  [2  [l  [0  C8  f|  *2  f3  ]  0  h  h  ]? 

[,  L  [,  [» o* f. 6 13  ]s  ]2  ]2  12  [,  [2  [,  [,c8f,f2f3  ]2  ]2  12  12  12 

Step  9: 

The  counter  has  reached  its  maximum  limit  since  2n+m-l  =  8  and  hence,  the  rule  applied 
would  be  [0  c  2„+m_,  ]  0  t ,  that  is  [0  c8  ]  0  -*■  t. 

Thus,  all  membrane  0  are  dissolved.  We  get, 

L  [3  [2 1, * t, t2 13  ]2  ]2  12  [3  h  [,ttit2f3  ]2  12  12 

[3  h  [, tt. fit3  ]?  ]2  ];  [3  L 1, * t, fi f,  ]?  ];  ]; 

[3  [2  [1  t  fl  t2  t3  ]  1  ]  2  ]  3  [3  [2  [l  t  fl  ]  1°  ]  2  ]  3 

[3  L  [.tfifata  ]?  ]“  ]3°  [3I2I.tf.f2f3  ]?  ]°2  ]?  ]2 

Step  10: 

Now,  the  rule 

[y  t,  ]  '■  — >  t| ,  if  x;  appears  in  clause  Cj,  1  <  i  <  n,  1  <  j  <  m,  and 

[ .  f,  ]  °j  —>  f, ,  if-  X|  appears  in  clause  Cj,  1  <  i  <  n,  1  <  j  <  m 

is  applied.  Clause  1  is  satisfied  by  all  except  fourth,  seventh  and  eighth.  Thus,  remaining 
membrane  1  are  dissolved.  We  get, 

[4  [3  [2 1  ti  t2 13  ]  °  ]  3  [3  [2 1  ti  t2  f3  ]  °  ]  3 
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[3  [2  t  tl  t3  ]°  ]°  [3  [2  [l  t1;l  f2  f3  ]  1  ]2  ]“ 


[3  [2  t  fi  t2  t3  ]  °  ]  3  [3  [2  t  fl  t2  f3  ]  °  ]  3 

[3  L  [1tflf2t3  ]?  ]°  ]3°  [3  [2  [x  t  f!  f2  f3  ]?  ]°2  ]°3 

Step  11: 

Clause  2  is  satisfied  by  third  and  fifth  truth  assignments,  so  the  corresponding  membrane 
2  are  dissolved  and  we  get, 

[4  [3  [2 1  fi  t2 13  ]  °  ]  3  [3  [2 1  ti  t2  f?  ]  2  ]  3 

[3 1  ti  f2  t3  ]°  [3  [2  [i  t  ti  f2  f3  ]  1  ]°2  ]° 

[3 1  fi  t2 13  ]  3  [3  [2 1  fi  t2  f3  ]  °  ]  3 

[3  [2  t fi f2 13  ]?  ]°  ]3°  [3  [2  [, t fl f2 f3  ]?  ]“  IS  IS 

Step  12: 

Clause  3  is  satisfied  by  only  third  truth  assignment  and  thus,  one  copy  of  t  is  left  free  in 
the  skin  membrane,  corresponding  to  the  truth  assignment  which  satisfy  the  formula. 

L  [3  [2 1  fi  t2 13  ]2  ]  3  [3  [2 1  ti  t2  f3  ]  °  ]°  tfi  f2t3  [3  [2  [,tti  f2f3  ]°  ]°2  ]° 

[3 1  fi  t2  t3  ]  3  [3  [2 1  fi  t2  f3  ]  °  ]  3  [3  [2  [1 1  fi  f2  t3  ]  1  ]  “  ]  3  [3  [2  [j  t  fi  f2 

f3  ]i°  ]°2  ]3°  ]S 

Step  13: 

t  will  be  sent  out  of  the  system  by  the  rule, 

[  t  1 0  — >  T  1+  t 

Lm+1  J  m+ 1  Lm+1  J  m+ 1 

Thus  the  formula  is  satisfied.  This  is  the  last  step  of  the  computation  because  no  further 
rule  can  be  applied.  The  skin  becomes  positively  charged. 
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Algorithm  to  solve  Hamiltonian  Path  Problem 

The  computational  steps  for  the  construction  of  the  P  system  for  the  Graph, 

G  =  ({1,2,3,4},  {{1,2},  {2,3},  {3,4}}) 

Here  n  =  4,  so  we  will  have  6  membranes,  labeled  from  0  to  5. 

Step  0: 

The  initial  configuration  is  given  as: 

LLLLLLcf  ];  p  p  p  p  p 

Step  1: 

The  counter  Co  is  incremented  using  Rule  1.  Membrane  0  is  divided  using  Rule  3. 

[5  [4  [3  [2  [1  [oCiai  ];  [0cia2  ]0  [0cia3  ]  +0  [0  Cia4  ]0  ] "  ]°  ]  3°  ]°  ]s 

Step  2: 

The  counter  Ci  is  incremented  to  c2.  Membrane  1  is  divided  according  to  Rule  9.  No  rules 
of  type  4  are  applied  to  ai  here  since  the  membrane  0  is  not  neutral. 

[s  [4  [l  [2  [|  [»  c2al  P  ]  :  [l  [»  C2a2  P  ] 

[,  [0 c2a3  ]S  p  [,  [0 c2a4  ]S  p  IS  IS  P  IS 

Step  3: 

The  counter  is  incremented  to  c3.  Rules  of  type  4  are  applied  to  a;.  Membrane  2  is  divided 
by  Rule  9. 

[5  [4  [3  [2  [1  [0  c3- 1  a2  ]  0  ]  1  ]  2  [2  [1  [0  c3-2.a3  ]  0  [0  c3-2.ai  ]0  ]  °  ]2 

[2  [1  [0  c3.3.a4  ]  o  [0  c3.3.a2  ]  0  ]  1  ]  2  [2  [1  [0 

CiA.ai  ]S  ]S  ];  IS  IS  IS 
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Step  4: 

The  counter  C3  is  replaced  by  C4.  Rules  of  type  4  are  again  applied  to  a2  and  a3  in  neutral 

membranes  0.  For  a in  non-neutral  membranes  0,  rule  of  type  4  cannot  be  applied; 

hence,  membrane  1  divides  instead.  Membrane  3  is  divided  using  the  rule  of  type  9. 

[5  [4  [.,[2  [,  [0  C4  .1.2.  a3  ];[,  c4  .1.2.  a,  p  ]f  ]S  p 

[3  [2  [,  [»  C4.2.a3  ]J  ];[,  [„  c4 .2.  a,  ]  °  ]  T  IS  P 

[3  [2  [,  [«  C4.3.04  ]“  ]r[,  [,  .3.  a2  ]  l  ]  ~  ]"  ]: 

[j  [2  [,  [,04.4.3.a4  ];[,  c4 .4.3.  a2  p  p  p  p  p  p 


Step  5: 

The  counter  C4  is  incremented  to  C5.  Membrane  4  is  divided  by  the  rule, 

r  r  i+  r  1  r  i+  i°  r  r  i°  i°  r  r  i°  i°  r  r  i°  t 

In  ln-\  An  — 1  L n  1  An-\ . \.n-\  An-\  An  In  ln-\  An-\  An  In  ln-\  An-\  An .  LnLw-1  An-\  An 


where  n  =  4. 

Membrane  1  is  divided  in  the  case  where  membrane  0  is  non-neutral.  Membrane  2  is 
divided  if  membrane  1  is  non-neutral.  Rules  of  type  4  are  applied  to  that  a;  which  is  there 
in  neutral  membranes  0. 

LLM  [,  [, c5. 1.2.83  is  p[,  [,  c5 .1.2. a,  ]s  ]r  p  is  p 

[4  [3  [2  [1  [0  05  .2.3.  84  ] ,  [,  C5  .2.3.  82  ]  0  ] ,  ]  2  [2  [1  [0  05  -2.1. 

a2  is  is  h  is  is 

[4  [3  [2  [.  [0C5  .3.4.  83  ]s  ] S  ]  2  [2  [,  [0 Cs  .3.2.  a3  ];U  .3.2. 

a,  ]S  ]?  ]S  IS  IS 

[4  [3  [2  [,  [,05.4.3.84  ]J  ];[,  [„  C5  .4.3.82  ]J  ]S  ]S  ]?  ]S  ]S 


Step  6: 

L  [4  [3  [2  [,  [.  c,  .1.2.3.84  ]  0  [0  C,.  1.2.3.82  ]S  IS  IS 

[2  [,  [,06.1.2.1.8,  IS  IS  Is  IS  IS 

L  [3  L  [,  [,  06-2.3.84  ]S  ]S[,  [„  06.2.3.82  ]J  ]f  ]S  P 


17 


[,  [,  [,  [o  C6.2-l-2.aj  ];[„(*  .2.1.2.  a,  p  P  P  p  P 

[,  [j  [2  [■  [o  C6.3A3.a4  ]J  [0  c6 .3.4.3.  a2  p  ]”  p  p 

[3  [2  [,  [0  c6  3.2.  a,  ]J  p[,  [0c63.2.a,  ]J  p  p  p  p 
l.l,L  [,  [«C6  4.3.4.  a3  ];  p  p 

[2  [,  [0  C6  4.3.2.  aj  ]  0  [0  C6  4.3.2.  a,  p  p  p  p  p  p 

Step  7 : 

[j  [4  [3  [2  [,  [.  C7.1.2.3.84  ]S  ]*[,  [„  c7. 1.2.3.  a2  p  p  p  p 
[3  [2  [,  [«  C,.  1.2.1.2.a3  ];  [,07.1.2.1.2.  a,  ],  ]?  ]!  p  ]] 

[4  [3  [2  [1  [.  C7.2.3.4.a,  ]  J  ]  ?  ]:[,[,[,  c,.2.3.2.a,  ] ;  [,  c,. 2.3.2. 

a,  ]„■  p  p  p  p 

[4  [3  [2 1,[. c7 2.1.2.33  ] ;  p[,  [,  c.2.1.2. a,  p  p  p  p  p 

[4 [3  [2 1,  [. c, .3.4.3.04  ] ;  p[,  [0c,.3.4.3.a2]s  p  p  ];  p 

[4  [3  [2  [l  [0  07.3.2.3.34  ],  [,  C7.3.2.3.32  ].  ],  p 
[2 [1  [0 c7.3.2.i.a2  ];  p  p  ];  ]; 

[4  [3  [2  [,  [.  C7  A3.4.3.  34  ]  ;  [,  c,  A3.4.3.  a2  ]„-  ]?  ]?  p 
[3  [2  [,  [,07  4.3.2.33  ]S  p[,  [,06.4.3.2.8,  ]J  ]-  ];  p  p  ]; 

Step  8: 

No  more  rules  are  possible  for  a;  in  membrane  0,  as  Ci  cannot  increase  further.  Thus,  Rule 
2  is  applied  and  the  counter  is  transformed  to  t  as  membrane  0  dissolves. 

L  L  [,[,[,  1.12.3.34  p  ];[,[,  1. 1.2.3. a,  p  p  p  ]; 

I,  [3  L  [,  t.i.2.1.2.3, .t.i.2.1.2.3,  p  ];  p  ]; 

[4  [3  [2  [,  12.3.4.3,  p  ];  p  [,[,[,  12.3.2.3,.  1.2.52.3-  p  p  p  ]J 

L  LU,  12-1-2.3,  p  p  [,[,  1.2.12.3,  ];  p  p  p 
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[.  L  L  [,1.3.4.3.34  p  ];[,[,  t  .3.4.3.  a2  p  p  p  p 
[4  [,[,[,  t.3.2.3.a4.t.3.2.3.a2  p  p  ]  J  [,  [2  [,  t.3.2.1.a2  P  P  p  ] 
[4  [3  [2  [,  t.4.3.4.3.  a4.t.4.3.4.3.  a,  p  p  p  ]J 
L  [3  [2  [,  1 4.3.2.  a3  ]”  p  [2  [,  t  .4.3.2.  a,  p  p  p  p  p 


Step  9: 

Rule  5  is  applied  and  all  the  afs  are  converted  into  i’s.  Membrane  1  dissolves  if  1  is 
present  according  to  rule  6.  Membrane  1  also  dissolves  if  ai  is  present  according  to  rule  7. 

[.,  [4  [,  [,  t.i. 2.3.4  p  p[3  [2 1 .1.2.3. 2  p  p  p 
[4  [3  [2 1 . 1 .2.1.2. 3. t .1.2.1. 2.1  p  p  ]; 

[4  [3  [2  [,  t.2.3.4.3  p  p  p  p 
[,  [3  [2 12.3.2.3.1.2.3.2.1  p  p  p 
[4  [3  [2 12.1.2.3  p  p  [,  [2 1 .2. 1.2.1  p  p  p 
[4  [3  [2  [,  t.3.4.3.4  ]?  ]:  p  [,[,[,  t.3.4.3.2  ]?  ]!  p  ]J 
[a  [3  [2  [,  t.3.2.3.4.t.3.2.3.2  p  p  p  p 
[4  [3  [2  13.2.1.2  ]?  p  ]J 
[4  [3  [2  [,  t.4.3.4.3.4.t.4.3.4.3.2  ]?  p  p  p 
P  [3  [2  [,14.3.2.3  p  ];  p  [,  [,  t  .4.3.2. 1  p  p  p  p 


Step  10: 

Membrane  2  dissolves  if  2  is  present. 

[5  Us  t-1-2-3-4  ]S  ]S  [4  [3  1. 1.2.3. 2  ]°  ]° 
[4  [3  t  .1.2. 1.2.  3.  t  .1.2. 1.2.1  ]°  ]° 

[4  [3  [2  [1  t.2.3.4.3  ]?]2  IS  12 

[4  [3  t.2.3.2.3.  t.2.3.2.1  ]  3°  ]  “ 
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[„[,  1 2. 1.2.3  p  p  [4  [3 1. 2. 1.2.1  p  ]; 

[4  [3  [2  [.  t.3.4.3.4  p  p  p  ];  [,[,[,[,  t. 3.4.32  p  ] 

[4  [,  [4  [,1.3-2.3  4.1.3.2.3.2  p  p  p  p 
[4  I  1-3  2.1.2  |  ] 


[4  [j  I.  [,1.4.3.4.3.4.1.4.3.43.2  p  p  ].  ] 

1. 1  r.  r.  14.3.2.3 1?  l"  1 


p  i 


Step  11: 

Membrane  3  dissolves  if  3  is  present. 

[,  [4  ti. 2.3.4  p  [4 1.1.2.3. 2  p  [4 1 .1.2.I.2. 3. t .1.2.1.2.1  p 

[4  [3  [4  [,  123.4.3  ]•  ];  ]•  ];  [4  12.3.2.3.1.2.3.2.1  ]J 
[4  1.2.123  |.  [4[  t  .2. 1.2.1  p  p 

[4  L  [4  [,  13.4.3.4  p  p  p  ];  [4[3[2[,  t.3.4.3.2  p  p  p  ]; 

[4  [3  L  [,13.2.3.4.1.3.2.3.2  p  p  p  p 

[4  13-2.1.2  ];  [4[3[2[,  1.4.3. 4.3.4.1.43. 4.3.2  p  |  |  p 

L  [3  L  [,  14.3.2.3  p  p  p  [I  [4  1.43.2.1  p  p 

Step  12: 

Membrane  4  dissolves  if  4  is  present. 

[5  t.1.2.3.4  [4  t.1.2.3.  2  ]°4  [4  1. 1.2. 1.2.  3. 1. 1.2. 1.2.1  ]4 

[4  [,  [4  [,  1.23.4.3  [MS  IN!  [4  1.2.3.23. 1.23.2.1  p 
[4  1.2.123  ];  [4[3  t  .2. 1.2.1  in; 

[4  [3  [4  [,1.3.43.4  p  p  p  ];  [4[3[2[,  13.4.3.2  p  p  p  p 

[4  [,  [4  [,  132.3.4.1.3.23.2  p  p  p  p 

[4  132.1.2  p  [4[3[2[,  1.4.3.4.3.4.1.43.4.32  p  p  p  p 
L  [3  [4  [,  14.3.2.3  p  p  p  p  1.43.2.1  p 
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Step  13: 

There  are  two  copies  of  t  in  the  skin  membrane.  One  of  the  copies  of  t  goes  out  of  the 
skin  membrane.  The  computation  stops,  because  there  are  no  more  rules  applied.  Thus, 
the  Hamiltonian  Path  Problem  is  solved. 
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